fxy=zeros(5,5);
fxy(2:4,2:4)=1;
subplot(3,2,1),imshow(fxy),xlabel('(b) 原始图像');
subplot(3,2,2),stem3([1:5],[1:5],fxy,'.'),axis([0 6 0 6 0 1]),xlabel('(c) 原始图像三维示意图');

Fuv=fft2(fxy);
ABSFuv=abs(Fuv);
subplot(3,2,3),imshow(ABSFuv,[]),xlabel('(f) 傅里叶幅度谱图像');
subplot(3,2,4),stem3([1:5],[1:5],ABSFuv,'.'),axis([0 6 0 6 0 9]),xlabel('(g) 傅里叶幅度谱三维示意图');

ShiftABSFuv=fftshift(ABSFuv);
subplot(3,2,5),imshow(ShiftABSFuv,[]),xlabel('(i) 对称傅里叶幅度谱图像');
subplot(3,2,6),stem3([1:5],[1:5],ShiftABSFuv,'.'),axis([0 6 0 6 0 9]),xlabel('(j) 对称傅里叶幅度谱三维示意图');